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Failure properties of loaded fiber bundles having a lower cutoff in fiber threshold 

distribution 

Srutarshi Pradharjj and Alex Hanseitj 
Department of Physics, Norwegian University of Science and Technology, N-7491 Trondheim, Norway 

Presence of lower cutoff in fiber threshold distribution may affect the failure properties of a bundle 
of fibers subjected to external load. We investigate this possibility -both in a equal load sharing 
(ELS) model and in local load sharing (LLS) one using analytic as well as numerical methods. In 
ELS model, the critical strength gets modified and beyond a certain lower cutoff level the whole 
bundle fails instantly (brittle failure) after the first fiber ruptures. Although the dynamic exponents 
for the order parameter, susceptibility and relaxation time remain unchanged, the avalanche size 
distribution shows a gradual deviation from the mean field power law. A similar 'instant failure' 
situation occurs in LLS model at a lower cutoff level which reduces to that of in the equivalent ELS 
model at higher (high enough) dimensions. Also, the system size variation of bundle's strength and 
the avalanche statistics show strong dependence on the lower cutoff level. 

PACS numbers: 



I. INTRODUCTION 

Critical behavior of the fracture-failure phenomena in 
disordered materials has attracted wide interest these 
days y). Among several model studies, the fiber bun- 
dle models (FBM) capture almost correctly the collective 
static and dynamics of fracture-failure in loaded materi- 
als. The two different versions of FBM have been stud- 
ied much. The equal load sharing (ELS) model 0, Q 
considers democratic (equal) sharing of applied load on 
the bundle, whereas in local load sharing (LLS) model 
0, |5j , only the nearest neighbors support the terminal 
load (stress) of a failed fiber. Experimentally it has 
been observed 0, 0, la| that disordered material under 
increasing load shows well defined power laws in terms 
of acoustic emissions prior to the global rupture. Such 
a power law in burst avalanches has been achieved ana- 
lytically (and verified through simulations) by Henrmcr 
and Hansen |!j,ll(| in ELS model. It is also known since 
several decades that the static ELS model has a critical 
point HHIll, E3, [H,[l5J], i.e., at a critical strength (a c ) 
the bundle shows a phase transition from a state of par- 
tial failure (for a < a c ) to the state of total failure (for 
a > a r )- This failure dynamics has been solved analyti- 
cally |16t fl7| which explores the critical behavior through 
the power law variation of order parameter, susceptibility 
and relaxation time. Also the mean-filed universality of 
ELS model has been established recently Qj| ■ The exten- 
sive studies on LLS model |J,Ia| suggest that the strength 
goe s to 'zero' value as the bundle size approaches infinity 
U3 l2fi , l2lj , -thus excludes the possibility of any criti- 
cal behavior. Another important observation is that no 
universal power law asymptotics exists for the avalanche 
statistics 10] in LLS model. Attempts have also been 
made to study both, the ELS and LLS models in a single 
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framework introducing adjustable load sharing parame- 
ter |22t l'23l l2-l| and a crossover from mean- field (ELS) to 
short-range (LLS) behavior has been reported. 

So far in the FBM studies, people mainly consid- 
ered different fiber threshold distributions starting from 
zero threshold. However, in reality every element (fiber) 
should have a finite (non-zero) strength threshold due 
to the cohesive force among the constituting molecules. 
Therefore the idea of a lower cutoff in fiber threshold 
distribution would be most welcome. Andersen et al |l3| 
considered first a lower cutoff in fiber threshold distri- 
bution and established the 'tricritical behavior' in the 
mean-field (ELS) fiber bundle model. Such distributions 
with lower cutoff have also been considered to study the 
nonlinear response in ELS mode [I3] and to establish the 
universal behavior [18( of ELS failure dynamics. The ex- 
clusion of weaker fibers not only enhances the ultimate 
strength of the bundle, it can affect the failure proper- 
ties of the bundle. To investigate such possibilities -in 
the present work we consider ELS and LLS models and 
proceed through analytic as well as numerical methods. 

We organize this report as follows: After a brief intro- 
duction (Section I) we study the effect of lower cutoff in 
ELS model (Section II) and in LLS model (Section III). 
The importance of such study and the physical signifi- 
cance of the observed results are discussed in the conclu- 
sion (Section IV). In the Appendix we apply our analytic 
formulations in two different situations of fiber threshold 
distribution. 



II. 'ELS' MODEL 

A. Solutions of the recursive dynamics for equal 
load increment 

We consider a fiber bundle model having N parallel fibers 
subjected to an external load or stress (load per fiber). 
The threshold strength of each fiber is determined by 
the stress value (<7th) it can bear, and beyond which it 



fails. We consider fiber threshold distribution to have a 
lower cutoff (<tl), i.e., a randomly distributed normalized 
density p(<J t h) has been chosen within the interval ul and 
1 such that 



p(a t h)d(Tth = 1- 



(1) 



We follow step- wise equal load increment [lfl [l7| till the 
total failure of the bundle. The breaking dynamics starts 
when an initial stress a (> ctl) is applied on the bundle. 
Fibers having strength less than a fail instantly reduc- 
ing the number of intact fibers and these fibers have to 
bear the applied load (ELS rule). Hence the effective 
stress (on intact fibers) increases and this compels some 
more fibers to break. These two sequential operations, 
the stress redistribution and further breaking of fibers 
continue till an equilibrium is reached, where either the 
surviving fibers are strong enough to bear the applied 
load or all fibers fail. 

The breaking dynamics can be represented by a recur- 
sion relation [l^ Il7l llSl] in discrete time steps: 



U t+1 = l-P(a/Ut);Uo = l; 



(2) 



where, Ut is the fraction of total fibers that survive after 
time step t and P(o~t) is the cumulative distribution of 
corresponding density p{o~th), 



P(.°t) 



p((Tth)d<Tth- 



(3) 



The time step indicates the number of stress redistribu- 
tion at a fixed applied load. 

At the equilibrium or steady state we get Ut+i — Ut = 
U* . This is a fixed point of the recursive dynamics and 
eqn. @ can be solved at the fixed point for some partic- 
ular strength distribution. 







Fig. 1: The uniform density of fiber strength having a lower 
cutoff (<Jl,) ■ 

We choose the uniform density of fiber strength thresh- 
old having a lower cutoff (Fig. 1) to solve the recur- 
sive failure dynamics of ELS model. Thus pi&th) has the 
form: 



p{vth) 



1 



1 - <7£ 



,0\L < Oth < 1. 



(4) 



The cumulative distribution becomes 

P(a t ) = H p(a th )da th = ^ ~ aL) 

Therefore Ut follows a simple recursion relation: 

1 



U t+1 = 



which has fixed points [171 Il8l ] : 
1 



1- 



u 



(5) 



(6) 



E/» 



2(1 -a L ) 



1± 1 



1/2' 



(7) 



where 



4(1 -a L ) 



(8) 



Beyond this critical strength (er c ) the whole bundle fails 
instantly. The solution with (+) sign is the stable one, 
whereas the one with (— ) sign gives unstable solution 

dm. 



(B) 




U 



U 



Fig. 2: The graphical solutions of eqn. JSJ : Straight lines 
represent the fixed points. In (A) ox = 0.3, cr c = 0.357, solid 
curve (a < <r c ), dashed line (a = cr c ), dotted line (cr > cr c ). In (B) 
ox = 0.5f , solid curve (cr = 0.52.), dashed line (<r = 0.55). 

It is obvious that the critical strength (ov) cannot be 
less than cxl- As a c is a critical point (see fla. \ll.\lcil ). 
there should be some critical exponents associated to a c . 
If none of the fibers fail (a < o~l), we cannot define order 
parameter, susceptibility, relaxation time etc. which show 
critical behavior of the failure dynamics. Therefore, at a c 
the bundle should be in a partially broken (stable) state. 
Putting the condition a c > o~l in eqn. (JSJ, we get the 
upper bound of lower cutoff: ctl < 1/2. We can verify 
that for ol > 1/2 the recursion (eqn. |BJ does not give a 
stable fixed point except U* — 0. Also, putting er^ > 1/2 
in fixed point solution (eqn. [JJ, we get U* > 1, which is 
unrealistic. Therefore the critical strength of ELS model 
is bounded by an upper limit: a c < 1/2. We present 
graphical solutions (Fig. 2) of the recursion relation (eqn. 
EJl for a L < 1/2 (A) and a h > 1/2 (B). Clearly, we cannot 
get a fixed point (U t +i = U t ) in (B). 



From the solution (eqn. UJ we can obtain the order 
parameter (O), susceptibility (%) and relaxation time (r) 
[HI 03 03 of the failure process: 

= U*(a)-U*(a c )^(a c -a)- a ;a=^ (9) 



X 



dU*{a) 



da 



(a c -a)^;(3=- 



1 



oY 



(10) 



(11) 



Therefore the variation of order parameter, susceptibility 
and relaxation time remain unaffected by the presence of 
lower cutoff. 



B. The 'instant failure' situation in weakest fiber 
breaking approach 

Now, we follow the weakest fiber breaking approach 
[3, [9] : The applied load is tuned in such a way that only 
the weakest fiber (among the intact fibers) will fail after 
each step of loading. We first find out the extreme condi- 
tion when the whole bundle fails instantly after the first 
fiber ruptures. As the strength thresholds of N fibers are 
uniformly distributed within <jl and 1, the weakest fiber 
fails at a stress ctl (for large N). After this single fiber 
failure, the load will be redistributed within intact fibers 
resulting a global stress <jf — N<tl/(N — 1). Now, the 
number of intact fibers having strength threshold below 
a f is 



ol 



NP(a f ) =N [ f p(a th )da th = ^^ 



(12) 



The stress redistribution can break at least another 
fiber if NP(<jf) > 1 and this 'second' failure will trigger 
another failure and so on. Thus the successive breaking of 
fibers cannot be stopped till the complete collapse of the 
bundle. Clearly, there cannot be any fixed point (critical 
point) for such 'instant failure' situation. Putting the 
value of Of we get 



(i-O 



>1; 



which gives 



o L > 



(N-l) 
(2N- 1)' 



(13) 



(14) 



For large TV limit, the above condition can be written as 
c_l > 1/2. Therefore, the condition to get a fixed point 
in the failure process is ox < 1/2. 
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Fig. 3: The total step of load increase (till final failure) is 
plotted against a^ for a ELS model having 50000 fibers. The dot- 
ted line represents the analytic form (eqn. 1171 . triangles are the 
simulated data for a strictly uniform strength distribution and the 
circles represent the data (averages are taken for 5000 samples) for 
a uniform on average distribution. 

We can also calculate how many steps are required to 
attain the final catastrophic failure for ctl < 1/2. Let 
us assume that we have to increase the external load x 
times before the final failure. At each step of such load 
increment only one fiber fails. Then after x step the fol- 
lowing condition should be fulfilled to have a catastrophic 
failure: 



r *i[l+l/(N-x)] 

N I p(<rth)d*th > 1. 



where 



0"l = 0"L 



x(l - a L ) 

N 



The solution gives 



x = ^(l '< 



2 V l-*£ 



(15) 



(16) 



(17) 



The above equation suggests that at <tl = 1/2, x = 0. 
But in reality we have to put the external load once 
to break the weakest fiber of the bundle. Therefore, 
x = 1 for <j l > 1/2 (Fig. 3). To check the validity 
of the above calculation we take 'strictly unifrom' and 
uniform on average distributions of fiber strength. In 
our 'strictly unifrom' distribution the strength of the k- 
th fiber (among N fibers) is ax + (1 — (TL)k/N. We can 
see in Fig. 3 that the 'strictly uniform distribution' ex- 
actly obeys the analytic formula Ijl7|l but the uniform 
on average distribution shows slight disagreement which 
comes from the fluctuation in the distribution function 
for finite system size. This fluctuation will disappear at 
the limit N — > oo where we expect perfect agreement. 



C. Avalanche size distribution 

During the failure process 'avalanches' of different size 
appear where simultanious failure of numbers of fiber 



is termed as 'avalanche'. To investigate whether the 
avalanche size distributions depend on the lower cutoff 
or not, we go for a numerical study. The result (Fig. 
4) demands that for small avalanche sizes the distribu- 
tions show a gradual deviation (depend on ox) from the 
mean- field result, although the big avalanches still follow 
the mean-field power law (exponent value —5/2) as an- 
alytically derived by Hemmer and Hansen ;9]. We have 
checked this result for several system sizes and the above 
feature remains invariant. We should mention that ac- 
cording to eqns. (37-39) of Ref. .9], the analytic treat- 
ment is valid for o^ < 1/2 as 1/2 is the maximum (upper 
limit of the integration in eqn. (37)) of the stress-strain 
curve. Presence of ox cuts out the lower part of the 
stress-strain curve where the small avalanches are most 
likely to happen. Therefore, smaller avalanches get re- 
duced in number with the increase of ox (Fig. 4). At 
the limit ox — > 1/2, the avalanche distributions seem to 
follow a new power law with exponent —3/2, which can 
be explained as: The fluctuation in threshold distribu- 
tion gives rise a scenario like the unbiased random walk 
of the bundle's strength around the maximum 1/2, which 
in turn results exponent —3/2 in avalanche distribution 
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Fig. 4: The avalanche size distributions for different values of 
ox: N = 50000 and averages are done over 10000 sample. We have 
drawn two power laws (dotted lines) as reference lines to compare 
our numerical results. 



neighbors (intact) having strength threshold below o/ 
is (nn)faii = nP{<jf) (see eqn. JUJl). Putting the value 
of P(<Jf) and <jf we finally get 



(nn)f a a 



(1-vlY 



(18) 



If (nn)f a ii > 1, then at least another fiber fails and this is 
likely to trigger a cascade of failure events resulting com- 
plete collapse of the bundle. Therefore, to avoid the 'in- 
stant failure' situation we must have (nri)f a u < 1, from 
which we get the upper bound of o^ : 



(JL < 



(19) 



As the above condition does not depend on the coordi- 
nation number n, at any dimension the whole bundle is 
likely to collapse at once for o^ > 1/2. It should be men- 
tioned that LLS model should behave almost like ELS 
model at the limit of infinite dimensions and therefore 
the identical bound (of ox) in both the cases is not sur- 
prising (see Appendix). We numerically confirm (Fig. 5) 
the above analytic argument (eqn. I19|) in one dimension. 
When average step value goes below 1.5 line, one step 
failure is the dominating mode then. We can find out 
the extreme limit of ox when all the nearest neighbors 
fail after the weakest fiber breakes. Then the LLS bundle 
collapses instantly for sure. Setting [nn)f a il = n we get 
the condition ox > n/(l+n), where stress level of all 
the nearest neighbors crosses the upper cutoff 1 of the 
strength distribution. Clearly such failure is very rapid 
(like a chain reaction) and does not depend on the shape 
of the strength distributions, except for the upper cutoff. 
Also as n increases (ELS limit) - ox for instant failure as- 
sumes the trivial value 1 . Similar sudden failure in FBM 
has been discussed by Moreno et al 25] in the context of 
a 'one sided load transfer' model which is different from 
the true LLS model -we consider here. 



III. 



'LLS' MODEL 



A. The 'instant failure' situation 

Now we consider LLS model with uniform fiber thresh- 
old distribution having a lower cutoff ox (Fig. 1). We 
shall present a probabilistic argument to determine the 
upper limit of ol , beyond which the whole bundle fails at 
once. Following the weakest fiber breaking approach the 
first fiber fails at an applied stress ox (for large N). As 
we are using periodic boundary conditions, the n near- 
est neighbors (n is the coordination number) bear the 
terminal stress of the failing fiber and their stress value 
rises to o/ = ox(l + 1/n). Now, number of nearest 
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Fig. 5: Numerical estimate of the upper bound of ox in LLS 
model: For ox > 0.5 the average step values go below 1.5, i.e., the 
bundle fails at one step in most of the realisations. 



B. Strength of the bundle 

The local load sharing (LLS) scheme introduces stress 
enhancement around the failed fiber, which accelerate 
damage evolution. Therefore, a few isolated cracks can 
drive the system toward complete failure through growth 
and coalescence. The LLS model shows zero strength (for 
fiber threshold distributions starting from zero value) at 
the limit N — > oo, following a logarithmic dependence on 
the system size (TV) [3 IM EH • Recently Mahesh et al 
26] have proposed a probabilistic method of finding the 
asymptotic strength of bundles in LLS mode. Now for 
threshold distributions having a lower cutoff (ox), the ul- 
timate strength of the bundle cannot be less than ul ■ For 
such a uniform distribution (Fig. 1), we perform numer- 
ical simulations to investigate the system size variation 
of bundle's strength. We observe (Fig. 6) that as o\l 
increases the quantity (strength-err, ) approaches zero fol- 
lowing straight lines with 1/N, but the slope gradually 
decreases -which suggests that the system size depen- 
dence of the strength gradually becomes weaker. 
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Fig. 6: The (strength — ox) is plotted against 1/JV for different 
ox values: 0.3 (square), 0.35 (circle), 0.4 (up triangle), 0.45 (down 
triangle), 0.5 (star). All the straight lines approach value as 

N — > oo. 



Avalanche size distribution 
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Fig. 7: Avalanche size distribution for several ol in LLS model: 
N = 20000, averages are taken over 10000 samples. Number of 
different avalanches decreases with the increase of ox value. 



Due to faster crack-growth, LLS model shows different 
avalanche statistics than that of in ELS model. The nu- 
merical study of Hansen and Hemmer [jj suggests an ap- 
parent power law having exponent —4.5 in the avalanche 
distribution. Later, Kloster et al [10J have shown ana- 
lytically that for flat (uniform) distribution LLS model 
does not have any power law asymptotics in avalanche 
statistics. We numerically study the avalanche distribu- 
tion in LLS model for different <tl values (Fig. 7). We 
observe a similar deviation (lowering) of the distribution 
function for the smaller avalanche sizes as in case of ELS 
model. Also, the number of different avalanches gets re- 
duced (tail of the distribution disappears) with the in- 
crease of &L . This occurs due to damage localization [27| 
which ensures faster collapse of the bundle. In Fig. 7 
we can see that for ol = 0.5, avalanches of size 1 is the 
only possibility before total failure and their count is al- 
ways less than 1, which clearly indicates the dominance 
of 'instant failure' situation. 



IV. CONCLUSION 

A lower cutoff in fiber threshold distribution excludes 
the presence of very weak fibers in a bundle. The weaker 
fibers mainly reduces the strength of a bundle. But 
in practical purpose we always try to build stronger 
and stronger materials (ropes, cables etc.) from the fi- 
brous elements. Therefore this situation (exclusion of 
weaker fibers) is very realistic. The failure dynamics 
of ELS model almost remains unchanged in the pres- 
ence of such lower cutoff (ox), whereas the avalanche 
size distributions show a systematic deviation (for small 
avalanches) from the mean field nature. At the limit- 
ing point (er L — » 1/2), we get a new power law (expo- 
nent —3/2) in avalanche distribution which can be ex- 
plained from random walk statistics [T^] ■ In LLS model 
the avalanche statistics show drastic change with the in- 
crease of a i,. In both the models, the lower cutoff be- 
comes bounded by an upper limit (ol < 1/2) beyond 
which the whole bundle fails at once, which has impor- 
tant consequence: It seems that the bundles show elastic- 
like response [l7| up to a^ = 1/2, above which they 
become perfectly brittle. We observe that the 'weak- 
est fiber breaking approach' y, an d the 'equal load 
increment approach' [16L Il7| give similar result (in ELS 
mode). In the equal load increment method, sometimes 
more than one fibers fail at the time of loading and this 
affects the whole failure dynamics, whereas the weakest 
fiber breaking approach ensures the single fiber (weakest 
among the intact fibers) failure at each step of loading. 
We consider the equal load increment method to be more 
practical from the experimental point of view. This ap- 
proach helps to construct the recursion relations [l(| [l7| 
which in turn show critical behavior |17L Hq of the fail- 
ure process. The 'instant failure' situation is not limited 
to uniform threshold distribution, rather it is common in 
any type of distributions (see Appendix). It seems that 



the 'instant failure' represents the binary states of the 
bundle: 'intact' (1 or high) and 'completely broken' (0 or 
low) and here the bundle behaves like a classical 'switch' 
in response to external load. 
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B. LLS model 

The condition for 'instant failure' for LLS model is: 

r <r L [l+l/n] 



p{cr th )da t h > 1. 



(23) 



where n is the coordination number or the number of 
nearest neighbors. This gives 



*i> 



i 



(24) 



Appendix 

Case I : Linearly increasing density of 
fiber strength 




Fig. 8: The linearly increasing density of fiber strength having 
a lower cutoff (ox) . 

We consider a bundle of fibers with linearly increasing 
density of strength having the normalised form: 

p{oth) = t^2, ox < oth < 1. (20) 

1 a L 

We want to find out the bound of ctl beyond which 'in- 
stant failure' occurs in both -ELS and LLS models. 

A. ELS model 

Following the weakest fiber breaking approach, the 
condition for 'instant failure' is: 

r <r L [l+l/(N-l)] 

N / p(ath)da t h > 1. (21) 

J CTL 



which gives 



ai> 



N-l 

3N-1 



(22) 



Therefore, the bound (beyond which one-instant failure 
will occur) of the lower cutoff comes to be ul < l/v3 
for large N limit. 



Now, as dimension of the system increases n goes towards 
infinity. Hence the above condition gives the bound of the 
lower cutoff as ox < l/v3, which is identical to that of 
in the equivalent ELS case. 

Case II : Linearly decreasing density of 
fiber strength 




Fig. 9: The linearly decreasing density of fiber strength having 
a lower cutoff (ox) . 

Next we consider a bundle of fibers having linearly de- 
creasing density of strength with the normalised form: 



2(1 -a th ) 

P(<Jth) = -r. To", OL < Oth < 1. 

(1-^l) 2 



(25) 



A. ELS model 

Following eqn. f!21|). the condition for 'instant failure' 



<JL> 



N - 1 



(26) 



3JV-1 

which sets the bound of the lower cutoff to ox < 1/3. 
B. LLS model 

Following eqn. ( 123(1 . the condition for 'instant failure 



is: 



2a L 



IX -°l) 



(l -,*-£)>!. (27) 



which reduces to 2ox/(l — ox) > 1 for n — ► oo and sets 
the bound of the lower cutoff to o^ < 1/3. Again this is 
identical to that of in the ELS case. 
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